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A MULTIPLE -IMPULSE FUNCTION FOR ORBITAL TRANSFER 
AND ITS DERIVATIVES 


By James F. Kibler and Lawrence H. Hoffman 
Langley Research Center 

SUMMARY 

A multiple -impulse function is represented as a sequence of single-impulse func- 
tions. The single-impulse transfer which yields the velocity change required to transfer 
from a specified initial orbit to a partially specified final orbit is developed. Analytic 
derivatives of the function are obtained for use in optimization techniques. A four- 
impulse transfer is outlined. The analytic derivatives may allow more efficient optimi- 
zation than numerical derivatives do. 

INTRODUCTION 

The calculation of the impulsive velocity change required to transfer between two 
general orbits is a frequently encountered problem in celestial mechanics. Orbital trim, 
planetary deboost, and orbital rendezvous problems are examples which may require one 
or more transfer maneuvers to achieve the desired objectives. In most cases, it is also 
desirable to minimize the fuel expended during the maneuver. A comprehensive survey 
of previous work in impulsive transfers has been made by Gobetz and Doll (ref. 1). In 
particular, the two-impulse transfer between inclined elliptical orbits has been studied 
intensively by McCue (ref. 2) and Lee (ref. 3). However, in some cases, it is desirable 
to satisfy constraints on the intermediate or final orbits and thus additional impulses may 
be required. The purpose of this paper is to develop a multiple -impulse function which 
describes the transfer between two general orbits. In addition, for minimization pur- 
poses, the partial derivatives of the function are taken with respect to the independent 
variables. 

The technique is to consider the multiple -impulse transfer to be a sequence of 
coupled single impulses; thereby, the orbit elements of the intermediate conics are made 
independent variables. The equations which describe the single -impulse transfer are 
found to be a function of only four variables on the second orbit; thus, the single-impulse 
maneuver consists of transferring from a given orbit to an orbit which is only partially 
specified. The multiple -impulse function can then be minimized over the space defined 
by four independent variables from each intermediate conic. 



The impulse function is developed for a single -impulse transfer. Analytic deriva- 
tives of the impulse function are presented and, as an example, the derivatives required 
to minimize a four-impulse transfer are given. 

SYMBOLS 

a '* semimajor axis, km 

e o eccentricity 

f true anomaly, deg 

i inclination, deg 

n number of impulses 

P unit radius vector lying in orbit plane and pointing toward periapsis 

Q unit radius vector lying in orbit plane and pointing toward 90° true anomaly 

R unit vector along line of intersection between two noncoplanar orbits 

r magnitude of radius at a given true anomaly, km 

r,h,n radial, horizontal, and normal components of velocity 

V velocity magnitude, km/sec 

AV magnitude of change in velocity required to transfer between two orbits, 

km/sec 

W unit vector normal to orbit plane and completing the PQW right-handed 

triad 

0 angle between two velocity vectors, deg 

y flight-path angle, deg 

6 angle between two orbit planes, deg 
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I u gravitational constant, kmVsec^ 

O right ascension of ascending node, deg 

<jo argument of periapsis, deg 

Subscripts: 

0,1 refers to the first or second orbit, respectively 

T total 

Bars over quantities denote vectors. 

ANALYSIS 

Derivation of the Impulse Function for a Single -Impulse Transfer 

The single-impulse transfer between two noncoplanar intersecting Keplerian orbits 
can be determined analytically if five elements are known on the first orbit and four 
elements are specified on the second orbit. It is arbitrarily chosen to specify the classi- 
cal elements (aQ,e0,i0,WQ,fio) on the initial orbit; and the second orbit is specified by 
the orientation angles (ii,wi,0;[) and one of the shape parameters (for example, aj). 

Of course, other choices are possible and would result in a different formulation. 

The first requirement is to find the line of intersection of the two orbits which is 
uniquely determined by ip, Oq, ij, and Oj. It is convenient to perform a coordinate 
transformation from the XYZ set of axes to the PQW system as shown in sketch (a) . 



Sketch (a) - Transformation of XYZ to PQW. 
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In the PQW system, P is a unit vector pointing in the direction of periapsis, Q is a 
unit vector lying in the plane of the orbit and pointing toward a true anomaly of 90°, and 
W is a unit vector, normal to the orbit plane, which completes the right-handed triad. 

By rotating through the angles £2, i, and u>, the resulting transformation (ref. 4) is 


* ^ 

P = (cos w cos £2 - sin a) sin £2 cos i) i + (cos oj sin £2 + sin w cos £2 cos i) j 
+ (sin w sin i)k 

Q = (-sin a) cos £2 - cos cu sin £2 cos i) i + (-sin w sin £2 + cos co cos £2 cos i) j ) 
+ (cos w sin i)k 

W = (sin £2 sin i) i + (-cos £2 sin i) j + (cos i)k 


( 1 ) 


where i, j, and k are the corresponding unit vectors in the XYZ coordinate system. 

The unit vector R, along the line of intersection of the two orbits, is determined 
by 


R = 


W Q XW[ 
W 0 XWj 


where 

Wq = (sin ig sin £2g) i + (-sin ig cos £2g) j + (cos ig)k 

(normal to first orbit) and 

Wj = (sin ij sin £2j) i + (-sin ij cos £2j) j + (cos ij)k 

(normal to second orbit) are known from equation (1). 

Now, 

I r - r 

i J k 


Wg XWj = 

sin ig 

sin £2g 

-sin ig 

COS £2g 

cos 


sin ij 

sin £2^ 

-sin ij 

cos £2j 

cos 


( 2 ) 


Expanding this expression results in 
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w o x Wi = 


(3) 


-sin ig cos cos fig + sin ij cos ig cos fij 
-sin ig cos ij_ sin fig + sin i^ cos ig sin fij 
sin ig sin ij sin Afi 

where Afi = fij - fig. The magnitude of Wg x Wj is given by 

|W 0 x Wj| = ^sin2ig cos^i^ cos^fig - 2 sin ig sin ij cos ig cos ij cos fig cbs fij 
+ sin2i^ cos^ig cos^fij + sin^ig cos^ij sin^fig 
-2 sin ig sin ij cos ig cos ij sin fig sin fij + sin2ij cos^ig sin^fij 
+ sin^ig sin^ij sin^Afi^^ 


Simplifying this expression results in 


IWgXWjl = [l - ^cos ig cos ij + sin ig sin ij cos Afi)^j 


1/2 


(4) 


Substituting equations (3) and (4) into equation (2) yields the expression for a unit vector 
in the XYZ system along the line of intersection between the two orbits: 


j/ - (cos ig cos q + sin ig sin q cos Afi)^j 


1/2 


-sin ig cos q cos fig + sin q cos ig cos fij 
-sin ig cos q sin fig + sin q cos ig sin fij 
sin ig sin q sin Afi 


(5) 


There are two points where an impulsive maneuver can take place, that is, in the direc- 
tion of either positive or negative R. Consider first the case of the transfer point lying 
in the positive R direction. The transfer -point true anomaly on each orbit must be 
determined. 

For the first orbit 
P 0 .R = 


R 


cos fg 


or 


cos fg = Pg • R 


(6) 
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Substituting the expression for Pg from equation (1), and equation (5) into equation (6), 
and simplifying yields 


cos fg = 


-cos o>g(sin ig cos ij - sin ij cos iQ cos Afi) + sin Wq sin ij sin Aft 


jl - (cos ig cos ij + sin ig sin ij cos Aft)^ 


1/2 


( 7 ) 


In order to resolve the ambiguity involved with the cosine, an expression for sin fg is 

.9 

required. 

By definition, 


P 


0 


X R = Wg 



R 


sin fQ 


or 

Wg sin fg = Pq x R 
From equation (2), 


Wq sin fg = Pg X 


W Q X W x 


WqXWj 


Pq X (Wp X W t ) 
Wq XWj 


t» 


Using the identity 

A x (b x c) = (A • C)B - (a • b)c 
where A, B, and C are any vectors yields 

,5, (Po ' - (h ■ *0)»1 

Wg sin fg = — 

|Wg X Wi 


But Pg • Wg = 0; therefore, 


Wg sin fg = 



0-*l) 

W ( 

3 


O 

X 



or 
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( 8 ) 


sin f n = tt 


Pq-Wj 


w 0 x w x 


Substituting equation (4) and the expressions for Pg and Wj from equation (1) into 
equation (8) and simplifying yields 


sin fg = 


sin u>g(sin iQ cos ij - sin ij cos iQ cos AO) + cos o>q sin ij sin AO 


l-( 


cos iQ cos ij + sin iQ sin ij cos 


AO)' 


nl/2 


O) 


In like manner, the equations defining the true anomaly of intersection on the second 
orbit are 

cos wRcos iQ sin ij - sin ig cos ij cos AOJ + sin u>i sin ig sin AO 
cos fj = — (10) 


1 - (cos ig cos i^ + sin ig sin ij cos Afi^ 


1/2 


and 


-sin o>i(cos ig sin i^ - sin ig cos ij cos Aft) + cos sin ig sin AO 

sin fj = — (11) 


I-( 


. r 

cos ig cos ii - sin ig sin ij cos AO)^ 


1/2 


The second point where an impulsive maneuver can take place is in the direction of nega- 
tive R. The true anomaly of intersection on the first orbit is given by 


cos fg = Pg • (-R) 


and 


^g sin fg = Pg X (-R) 

Substituting and simplifying as in equations (7) and (9) yields 

cos o>g(cos ij sin ig - sin ij cos ig cos AO) - sin wg sin ij sin AO 


cos fg = 


1 - (cos ig cos ij + sin ig sin ij cos Afij^ 


jl 1 / 2 


( 12 ) 


and 
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(13) 


sin fg = 


-sin o>g(cos ij sin ig - sin ij cos ig cos Aft) - cos u>g sin ij sin Aft 


1 - (cos ig cos ij + sin ig sin i^ cos Aft) 2 ] 


1/2 


Similarly, the equations defining the true anomaly of intersection on the second orbit are 


-cos ojj(cos ig sin ij - sin ig cos ij cos Aft) - sin sin ig sin Aft 
v cos fj = (14) 


i-0 


cos ig cos ij + sin ig sin ij cos Aft)* 


ill/2 


and 


sin fj = 


sin o>i(cos ig sin ij - sin ig cos ij cos Aft) - cos sin ig sin Aft 


1 - (cos ig cos ij + sin ig sin ij cos Aft) 2 


1/2 


(15) 


In other words, 

( f 0)-R ’ ( f 0) t R + 180 ° 

and 


('O-R * ( f 0 + R + l *°° 

For an impulsive maneuver to take place at the point of intersection, the magnitudes of 
the radii on the first and second orbit must be equal. Therefore, a constraint equation at 
the point of intersection is 

a 0 (l - e 0 2 ) ai (l - ei 2 ) 

1 + eg COS fg 1 + ej COS fj ' 

The elements ag, eg, fg, fj, and a^ or e^ are known. Therefore, equation (16) 
relates the unknown element in the second orbit ^ej or a^ to known parameters. 

Next, an expression for the magnitude of the vector difference in velocity (AV) 
between the two orbits is required. From sketch (b), using the law of cosines, 

(AV) 2 = v 0 2 + Vj 2 - 2VgV! cos 0 (17) 


where 0 is the angle between the two velocity vectors. 
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At the maneuver point, by using the equations 

11/2 


V = 




and 


r _ a (* - e2 ) 

1 + e cos f 

the expressions for the magnitude of the velocity on the first and second orbit are 

] 1/2 

i pi> ii t titn v/uo in t cn j 

v o = 



^1 + 2 eg COS fg 

+ e o 2 ) 


a 0 

1 - *o 2 ) 



(18) 


and 


Vf = 


M 

^1 + 2e 

L COS fj 

+ el 2 ) 


a l 

(> - e. 2 ] 

_ 


1/2 


(19) 


From sketch (c), by using the law of cosines for spherical triangles, the angle 
between the orbit planes 0 may be found from 

cos 6 = -cos iQ cos(l80° - ijJ + sin iQ sin(l80° - ijJ cos Afi 
or 

cos 8 = cos ig cos ij + sin ig sin ij cos A SI (20) 
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From sketch (d), the velocity vectors may be written as 


and 


V 0 



sin y Q 
cos y 0 


0 


Vl 


= V 1 


sin yj 

cos yj cos 6 
cos y-jL sin 9 


where Yq and y^ are the flight-path angles on the first and second orbits. The angle 
between the velocity vectors 0 may now be found from 


cos /3 = 


Vo ■ Vl 


|V 0 Vl| 


= cos 9 cos Yq cos y^ + sin y^ 


sin 


( 21 ) 
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Sketch (d) - Velocity vectors in the radial, horizontal, 
and normal coordinate system. 

Substituting the flight-path angle relations from reference 4 


sin y = 


e sin f 


(l + 2e cos f + e 2 ) 


1/2 


and 


cos y = 


1 + e cos f 


1/2 


(l + 2e cos f + e 2 ) 
and equation (20) into equation (21) yields 

(l + eg cos JoXl + cos f]j(cos i.Q cos + sin ip sin cos Aft) + (eg sin fg)(e^ sin f^) 


cos (3 : 


(l + 2eg cos fg + eg 2 )(l + 2ej cos fj + ^ 


(22) 


Substituting equations (18), (19), and (22) into equation (17) results in the required 
expression for AV : 


AV 


{— 77 ^ (* + 2e 0 cos f 0 + e 0 2 ) + —7 7 57 ( J + 2e l cos f l + e l 2 ) 

M 1 ' e o ) a i ( 1 - e i ) 

(l + eg COS fg)(l sin ig sin ij cos Aft) 


[ag(i - e 0 2 ) ai (l - ei 2 ] 1/2 
' 1/2 


+ ege^ sin f ( 


0 sin f l 


(23) 
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The set of equations (7), (9) to (16), and (23) are those required for an analytic solution 
to the problem of a single -impulse transfer between two noncoplanar orbits. Two solu- 
tions result from the positive and negative cases for R. In addition, if ej is a function 
of a^ in the constraint equation (16), two solutions arise from the quadratic in ej. All 
the multiple solutions may not be physically possible. For example, an arbitrarily spec- 
ified ai may result in a maximum radius on the second orbit which is less than the 
smallest radius on the initial orbit. That is, the constraint equation (16) cannot be sat- 
isfied. In addition, if the central body has a finite radius, a region of mathematically 
possible orbits must be rejected since they intercept the surface of the planet. In the 
event no practical solutions for the transfer exist, a more realistic set of variables 
(ai,ii,wi,fij} on the second orbit must be chosen. In the event of multiple practical 
solutions for the transfer, a choice may be based on mission considerations, such as 
choosing the transfer which results in the lowest AV. 


Partial Derivatives of the Single-Impulse Function 

In orbital transfer problems, it is usually desirable to minimize AV with respect 
to the control parameters. Many of the minimization algorithms require partial deriva- 
tives of AV with respect to the independent variables. Finite-difference approxima- 
tions to the partial derivatives are commonly employed. However, in complex multiple- 
impulse problems, it may take relatively large amounts of computer time to develop the 
approximations. If exact partial derivatives are available, the minimization algorithms 
are generally more efficient and less computer time is used. Therefore, analytic partial 
derivatives of equation (23) have been developed. It is assumed that a, e, i, w, and 
£2 are known on the initial orbit and that a, i, o>, and ft are known on the final orbit. 
Eccentricity of the final orbit is a dependent variable defined in the constraint 
equation (16). 

For convenience, define the following intermediate quantities as 


1-1 


A 1 = M a l(l - ®1 2 ) 


-1 


Bj = 1 + 2ej cos fj + ej/ 


AO = h ^ q(i - e 0 2 ) ' 

Bq = 1 + 2eg COS fg + eg 2 
C = -2/i[ag(l - e 0 2 ) ai (l - ei 2 | 1/2 

D = (l + eg COS fg X 1 + e i cos fj^cos ig cos ij + sin ig sin ij cos Aft^ 
+ egej sin fg sin f j 
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Thus, equation (23) can be written as 

AV = (AqBq + + CD) 1//2 


or 


(AV) 2 = AgBg + + CD 

The required derivatives can be expressed in the form 


3AV 1 3(AV) 2 

3 a 2AV da 


2AV 


b„^ + a 0 3B ° 


da 


da 


+ Bi 


3Aj 

da 


3Bi 

A 1 —± + 

da 


3D 

3 ( 3 ! 


+ D^ 

da) 


(24) 


where a represents any one of the variables ag, eg, ig, o>g, fig, a^, ij, Wj, 
and fij. 


It is necessary now to develop the partial derivatives of the following six interme- 
diate parameters: 


(1) 


3Ag 

3 Of ‘ 


Let 


Ag = /r|ag(l - e Q 2 ) 


n-i 


then 


3Ag 

9ag 



Ao 

*0 


Similarly, 


9Ap _ 2Ageg 
3e 0 “ 1 - e 0 2 

and 

9 Ap _ SAg _ 3Ag _ 3Ag _ 

3ig 3u)g 3 fig Sq^ 

where represents aj, i^, wj, and fij. 
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(2) 


9Aj 

8a 


Let 


A 1 = ^[ai(l - e^) 


1-1 


then 


= 2 (I e 1 a,-»(l - e,2)' 2 ifl = 2 ^ Al if 1 
9a 0 1 1 v 1 ' 9a 0 1 _ ei 2 9a 0 


Similarly, 


9Aj 2ejAj 9ej 
9e 0 1 - e^ 9 e 0 


9Aj 2e^Aj 9ej 

^Tl-e^ 

9Aj 2ejAj 9ej Aj 
9a l 1 - ei 2 **1 


9Aj 2ejAj 3ej 
aw l 1 - ei 2 9w l 


The 


9ej 

9q» 


remain to be determined. 


(3) 


9q» 


Let 


Bq = 1 + 2eg cos fg + e Q 2 


Then 


since 


9 Bp 

9a 0 


= 2e 0 


9 cos fg 

— = 0 

9ag 


9 COS fg 

3a 0 


= 0 


9Ai 2ejAj 9ej 


9ig 

1 - ei 2 

9ig 

9Aj 

. 2e l A l 

9ej 

3£2q 

1 - ej 2 

3S2q 

9A X 

2CJAJ 

9ej 

9il 

1 - ej 2 

9i l 

9A X 

2CJAJ 

9ej 


9^1 1 _ ej 2 9S2i 
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by inspection of equation (7). Similarly, 



+ cos fg) 


9Bq 

Sin 


= 2 e r 


9 COS fo 


31, 


0 


8B 0 9. 

awg * 2e 0 


9 COS fg 

9wg 


9Bq o 9 COS fg 
9f2 0 6 ° 9%) 



9Bq 

9i l 


= 2e 0 


9 cos f o 
ai l 


9Bg 

9ci>i 


= 0 


9Bg 

Sflj 


= 2e 0 


9 COS fg 
9f2j 


9 COS fg 

The — x remain to be determined. 

9q! 


9Bi 

(4) Let 


then 


f 

Bj = 1 + 2ej cos fj + e ^ 


9Bi 9 cos fi . 9ei 9ei 

- 2e i -s^r- - 2 ( cos f i + e i) - 2 ( cos f i + e i) ^ 


Bag 


9ag 


since 


9 cos fi 


”3ai 


0 


= 0 


by inspection of equation (10) . Similarly, 


9Bi 


9 eg ~ 2 ( Cosf l + e l)9iJ 


9Bi 


9 COS f 


” + 2 (cos 


9ig 


Sir 


9Bi . . oei o£>i o cos ii . . oei 

J- = 2(cos fi + e x )— MS = 2el SaT * 2 ( cos f l + e l) a ST 


9ej 


SB 


9 cos fi 


9ej 
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3Bi , . 3ei 

— = 2(cos fi + ej — 


3Bi 

— - = 2ei 

d(x) j 


3 cos 
3o)j 


9 ei 


+ 2(cos ti + t-~ 


3co^ 

3 COS fi 

The — - remain to be determined. 

3 O' 


9B 1 ■„ 3 COS f j 3 ei 

= 2e l— al^— + 2 ( cos f i + e l) 


3ij 

3B! 

3fii 


3 COS fi 

= 2ei an. 


, . 9 e. 

+ 2(cos fl + ei)^ 


(5) f : Let 


then 


Similarly, 


c = 

~2m|^o(i 

-e 0 2 ) 

3C 

9a o 

= M 

fol 1 - 

-e 0 2 )a 


- cl 

' e l 

_!!l 


-C\ 

/ " e i 

2 3a 0 

3C 

- d 

/ e x 

3ej 

9e o 


4 - e 

l 2 9e 0 

3C 

Cej 

3e^ 

9i 0 

1 - 

■ ej2 

3i 0 

3C 

Cej 

3ej 

3co 0 

1 

- e^2 

3w 0 

3C 

Ce l 

3ej 

3J2q 

1 

- ei 2 

30q 

3C 

-cf 

e l 

3ej 

3aj 

[ 

1 - e 

l 2 3 ai 

3C . 

Ce l 

3e^ 

9i l " 

" 1 - 

■i* 



1 - 1/2 


■3/2 


3ej 

3a0 
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9C 


9 (t) 


Cej 9ej 
2 9o)i 


1 1 ' e l 


3C _ Ce l 9e l 
9^1 \ - 0^2 


(6) |g: Let 


D = (l + eg cos f 0 )(l + cos f^cos ig cos ij + sin ig sin ij cos Aft;) 
+ egej sin fg sin fj 


Define the intermediate terms as 


Dn = 1 + 6n cos ffl 


Dj = 1 + ej cos fj 

D2 = cos ig cos ij + sin ig sin ij cos Aft 


then 


^ = (D 0 D 2 cos fj + e Q sin f 0 sin fj) ^ 

9D n L 9e l . r. 8e l 


— = D 2 |D 0 cos fj — + Di cos f 0 j + sin f 0 sin i x \e Q — + e x 


9D n 


Dj^cos ig sin i x cos Aft - sin ig cos i x ) + D 2 fe x 


9 cos f x 
9in 


9e 


+ cos fi — — 
1 9l r 


+ D 2 D l e 0 


9 COS f r 


9ir 


e 0 


/ 9 sin fi 9 sin fg\ 9e x 

e l( sin f 0 ai Q + Sin f l — aiT - ) + sin f 0 sin f l 


9ig 


9ig 


= °2 

9wq 


9 COS 


9ei 


9 cos f n 


D 0l 6 l-9^- +COS f l 9^) +D l e 0 


9o)fi 


+ e f 


9 sin f x 9 sin fg 


e i( sinf o-^^ + sin f l 


0 


9ei 

. , + sin fg sin fi -r — 

9d)g J u I 
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9D 

9 fig 


= D r 


/ 9 COS fl 

sin ig sin ij sin A fi + D2lei — — + cos fi 


9ei 


+ D l D 2 e 0 


9 COS fr 


9 fir 


e 0 


ei sin f 0 


9 sin f i 
9 fin 


+ sin fi 


9fi 0 

9 sin fn\ 9ei 

— — - — - + sin fn sin fi — — 
9fi o / 0 1 9^0 


^•7 = (d 0 D 2 cos fi + e 0 sin f 0 sin fi)^ 


9a i 


— = d 0 

9il 


. . / 9 cos fi 9ei 

D^sin ig cos ii cos Afi - cos ig sin iij + D2 (ei — — + cos fi -^ 7 — 


9 cos fn 


9 sin fi 


9 sin f, 


+ D i D 2 e 0 ^ + e o e il sin f o ~ aTj— + sin f l - 8 i 1 -j + sin f o sin f i 9 ^- 


9e, 


9D 

9o>i 


,D 2 


9 COS fl 3e l\ 8 cos 

D ol e l 8ui' + C0S fl 9^7 + Dl6 ° 8Wi' ' 


+ eg 


9 sin fi 9 sin fg\ 9ei 

e l( sin f 0 + sin f l + Sin f 0 sin f l 9717 


9D 

9fii 


= Dg 


-Di sin ig sin ii sin Afi + D 2 ei 


9 cos fi 9ei 

9fil + C0S fl 9fi7 


9 COS fg 

+ DiD 2 e 0 -gj^— + eg 


e x sin fg 


9 sin fi 
9fii 


+ sin fi 


9 sin f, 


9ei 

, + sin fn sin fi — — 

9«l J VI afij 


9 sin fn 9 sin ^ 

The and — are to be determined. 

da 9 q/ 

The preceding equations involve derivatives of eccentricity on the second orbit with 
respect to the elements on the first and second orbit. Rewriting the constraint equa- 
tion (16) yields 

a l(l - ei 2 )(l + e 0 cos f 0 ) = a 0 (l 
When E is defined as 

E = ag(l - eg 2 ) cos fj + 2aiei(l 

this equation can be differentiated implicitly to obtain the following derivatives: 


- eg 2 ^(l + ej cos i^j 
+ e Q cos f 0 ) 
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9e l 1 


3a r 


= |j(l - e 0 2 )(l + e l cos fj) 


= f a l(! " e l 2 ) cos f O + 2a O e o( 1 + e l cos f l) 


8 ei _ 1 

3i 0 ' E 

9e l _ 1 

9o>q E 

9e l _ 1 

3f2g E 


/ ON 3 COS fn / o\ 3 COS f 1 

~ e i - W 1 - e o 


- e l 2 ) 


3 COS f. 


3 <t)r 


/ 9 \ 9 COS fQ / 

a i e o(‘ - e i 2 )~5nJ - “oni 1 - V) 


3 COS f] 
3 fig 
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Finally, it is necessary to obtain derivatives of equations (7), (9), (10), and (11). (This 
method assumes that the transfer takes place in the direction of positive R. For nega- 
tive R (eqs. (12) to (15)) simply interchange the initial sign on the derivatives.) 

Define 

F = 1 - (cos ig cos ij + sin iQ sin ij cos 


By inspection of equation (7), 


3 cos fg 
9a 0 


3 COS fg 3 COS fg 
3eg 3aj 


3 COS fg 
3u)j 


3 COS fg 

8 H) 


-i^sin ij sin Afl(cos ig cos ij + sin ig sin ij cos AO) sin fgj 
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3 COS fg 


dC0 Q 

= sin fg 

3 COS fg 

ir. ■ ( 

3ftg 

~ FL in "IV 

3 COS fg 

1/ . . 

9i l 

= p( Sin ^ £ 


= sin ij (cos ig sin ij - sin iQ cos ij cos Aft) sin fg 


3 cos fn 3 cos f r 


30j 3ftg 

By inspection of equation (9), 


3 sin fg 3 sin fg 3 sin fg 3 sin fg 


Sag 


3e 


0 


3a i 


3uh 


= 0 


3 sin fg i r , .. 

— — = p|sin ij sin A ft [cos ig cos ij + sin ig sin ij cos A ftj 


= -cos f r 


3 sin fg 
9 ^0 

3 sin fg j 

3f2g F 


sin ij (sin ij cos ig - cos i^ sin ig cos Aft) cos fg 


3 sin fg i . 

— = - —(sin ig sin Aft cos fgj 

3 sin fg 3 sin fg 

3ftj 3ftg 


By inspection of equation (10), 
3 cos fi 


Sag 


3ig 


3ftg 


3 COS f 1 

3 COS f J 

3e 0 

II 

o 

3 

CO 

-?( sin *1 

sin Aft sin 

[ sin ‘o 

(sin ig cos : 


3aj 


= 0 


cos 


20 



3 COS fj 1 
9 COS fj 


= sin ig sin A£2 (cos ig cos ij + sin ig sin ij cos A sin f^ 


■= sin f i 


3 COS f 1 3 COS fj 


3fi 


1 


3 Or 


By inspection of equation (11), 

3 sin fj 3 sin f j 3 sin fj 3 sin fj 


9a r 


9e 


0 


dC0 r 


3a i 


= 0 


3 sin fi i . v 

— - = p( sin 11 Sin Afi C0S f l) 


3 sin fi i r , \ ~\ 

— — — = p[sin ig (sin ig cos ij - cos ig sin ij cos A Cl) cos fjJ 

3 sin fi I r / \ ~i 

— aT" — = “ plf in H) sin Af2(cos ^0 cos *1 + s ^ n k) sin H cos Afij cos f lJ 

3 sin f j 


8»f- " - C0S f l 


3 sin f j 3 sin fj 

30j 3J2g 


All these relations can be substituted into equation (24) to obtain the desired derivatives 
of AV. It is prohibitive to write these equations explicitly. However, they are readily 
programed for the computer. The analytic derivatives have been checked with good 
agreement against central -difference numerical derivatives. 


Generalization to Multiple -Impulse Transfer 

The multiple-impulse transfer is treated as a sequence of single-impulse transfers. 
The first orbit is specified by a, e, i, o>, and Cl, and the final orbit is specified by 
a, i, ci>, and Cl, eccentricity being a dependent variable. On each intermediate conic 
there are four control variables (a, i, w, and 0) and one dependent parameter (e). By 
choosing values of the control variables, the total velocity change is the sum of the single 
impulses. For minimization purposes, the derivatives of total AV with respect to the 
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control variables are obtained by chaining together the derivatives obtained in the pre- 
vious section for the single-impulse transfer. 

As an example, consider a four -impulse transfer problem. Let represent any 
of the four independent parameters of the ith intermediate orbit (i = 1,2,3). The total 
velocity change is the sum of the four impulses and may be written as 


AV'p = AVi + AV 2 + AV g + AV 4 


(25) 


where 

AVj = AVi^) 

AV 2 = AV 2 (o! 1 ,Q! 2 ,e 1 ) 
AV 3 = AV 3 ( 0 ! 2 ,a; 3 ,e 2 ) 
AV 4 = AV 4 (a 3 ,e 3 j 

and 


e l = e lK) 
e 2 = e 2( a l> a 2» e l) 
e 3 = e 3(p i 2> a 3> e 2) 

The partial derivatives of AV X with respect to the independent elements are 

9AV 3 /'9e 2 9e 2 9ej 

9e 3 l 9 ®! 9ej 9a?j 


9AV 4 9e 3 ^9e 2 9e 2 9e- 

+ 9e 3 9e 2 + 9ej 9a 

9AV t 9AV 2 /9AV 3 9AY 3 9e 2 ^ 
9a 2 9a 2 + \ 9a 2 + 9e 2 9a 2 y + 

9AV X 9AV 3 1 9AV 4 9AV 4 9e 3 
9a 3 9a 3 l 9a 3 9e 3 9 q! 3 



9AV 4 / 9e 3 


9e 3 9e 2' 


9e 3 \9 q! 2 9e 2 9a 2 y 


(27) 



9AV X _ 9AVj /9AV 2 9AV 2 9ej N 
9aj 9aj l da\ 9ei da\ 
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Thus, to relate any two consecutive orbits, the following derivatives are necessary: 

9AV|,i 9AV;,i 9AV;,i 9ej.i 9 e4.i 9e; 

— 5 , — 5 , — — - — , — T — , — — , and The derivatives are obtained from 

9ei 9 o?i ’ 9 a?i+i 9 ei ’ 9^ ’ 9^ 

equation (24) and AV i+ j represents the velocity change required to transfer from the 

ith to the (i + l)th orbit. The extension of this formulation to n-impulses simply involves 

n - 1 equations similar to equations (25) to (28). 

CONCLUDING REMARKS 

A multiple -impulse function has been developed as a series of single impulses. 
Analytic derivatives of the function with respect to the independent variables have been 
presented. An example of a four-impulse transfer is given and could be expanded to 
n-impulses by the same techniques. The development of analytic derivatives of this 
impulse function removes the dependency of optimization algorithms on numerical deriva- 
tives. Analytic derivatives may permit much more efficient optimization techniques. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., January 18, 1972. 
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